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Abstract 

The next-to- leading order three-jet cross section in hadron collisions is calcu- 
lated in the simplified case when the matrix elements of all QCD subprocesses 
are approximated by the pure gluon matrix element. The longitudinally- 
invariant k± jet-clustering algorithm is used. The important property of 
reduced renormalization and factorization scale dependence of the next-to- 
leading order physical cross section as compared to the Born cross section is 
demonstrated. 
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The structure of hadronic final states in high energy collisions may be described in 
terms of jet characteristics. Nowadays, jet cross section data are used both for precise 
quantitative tests of QCD — such as measurement of the strong coupling a s and QCD 
scale Aj^g — and looking for signs of new physics beyond the standard model. Both of 
these objectives have been studied extensively in e + e~ annihilation that is characterized 
by a fixed high energy scale, namely the machine energy. One would also like to analyze 
the data of scattering processes at the highest available energy scale, where the signs of 
new physics is expected to be the most profound. The highest scale is in general found in 
hadron collisions. Unfortunately, there are important ambiguities which limit our ability to 
perform high precision quantitative studies with the jet cross sections observed in hadron 
collisions. In experiments ambiguities arise from the question of how to define a jet and from 
the systematic uncertainties of jet energy measurements. On the theoretical side, apart 
from the ubiquitous and ever decreasing uncertainty in the parton density functions [0, 
there still remains uncertainty in the choice of renormalization and factorization scales, the 
magnitude of the higher order corrections and the question of how to match the theoretical 
and experimental jet definitions. The theoretical ambiguity coming from these points can 
be decreased if the next-to-leading order corrections are calculated. 

The most easily calculated next-to-leading order corrections to cross sections in hadron 
collisions are those to the inclusive one-jet and two-jet cross sections which have been avail- 
able for some time These results have already proven to be extremely important: the 
significant enhancement in the experimental one-jet cross section over the result of the the- 
ory in the range p!p > 200 GeV may be interpreted as signal of physics beyond the Standard 
Model §. 

In order to be able to calculate next-to-leading order corrections for more complex fi- 
nal states than the ones mentioned above — such as the three-jet cross section in hadron 
collisions — , two obstacles had to be overcome. Firstly, there was the issue of loop ma- 
trix elements which only recently have become available for all five-parton processes |§-|3| 
necessary for a three-jet analysis. Secondly, the algorithm for the cancelation of infrared 
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divergences applied in ref. [|J were not directly applicable in processes with more complex 
final state kinematics. Recently, a number of general schemes have been proposed for achiev- 
ing the cancelation of final state infrared singularities and mass factorization both in the 



framework of "phase space slicing" method |9| and "subtraction" method flT0|-|T2|. 

In this letter we present a brief summary of an analysis of three-jet cross sections using 
the algorithm of ref. [12|, but in the simplified case when all squared matrix elements are 



approximated with that of the pure gluon subprocess. Thus the results shown are intended 
only for demonstrating the applicability of the subtraction scheme of ref. [12| in the case 
of hadron collisions rather than a serious theoretical description of the data. We anticipate 
however, similar conclusions as those drawn here will hold once the complete analysis with 
quarks included is finished. 

According to the factorization theorem, the next-to-leading order infrared safe physical 
cross section at order a[ N+1 ^ is a sum of two integrals, 

<y = X[2 — > iV] + X[2 — > N -\- 1], (1) 

where in the case of the pure gluon approximation for the squared matrix element, these 
integrals have the form 

J[2 -> n]= J dx A fcs(g,x A ) J dx B f e s{g,x B ) (2) 



x^^^ / dT^(p1, ..., p %)(\M(g + g^ ng)\ 2 )SM, . . . ,p£). 
zxaXbs n\ J 

In this equation dr( n ) is the usual n-particle phase space measure. There are two possibilities 
for the choice of the effective gluon densities. One can either imagine colliding glueballs || , 
use a reference gluon density at a fixed /io scale and evolve it to other scales with Nf = 0, 
or alternatively, one can use the effective gluon density 
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fes{g,x) = f{g,x) + ~J2if(Vi x ) + f(^ x )]- (3) 
9 ? 

Our choice will be the latter one. There is still the question of what sort of strong coupling 
a s (fi) one should use in the pure gluon (\M.(g + g — > ng)\ 2 ) squared matrix elements. 



We use the two-loop formula for the strong coupling with Nf = and the QCD scale 
parameter in the modified minimal subtraction scheme, A^ chosen to be 1600 MeV so as 
a s (50 GeV) ~ 0.13 is consistent with the value of a s (50 GeV) observed in the real world with 
quarks. In this way we ensure that the relative size of the next-to-leading order correction is 
similar to that in the full QCD case. Note however, when we compare the order a[ N+1 ^ cross 
section to the results of the Born-level calculation, we compute the Born cross section using 
the one-loop formula for a s with the = 1100 MeV, which makes a s (50 GeV) about 15 % 
and the Born-level three-jet cross section about 50 % bigger. 

The function S n in eq. (0) represents the physical quantity to be calculated. Among the 
numerous infrared safe physical quantities one can calculate with the technique presented 



ref. [12], an explicit example, that we use in the present analysis, is the next-to-leading order 



three-jet cross section in hadron collisions defined using the longitudinally-invariant k± jet- 



clustering algorithm [TJ]]. For hadron collisions the jet-clustering algorithm is a two-stage 
process, each characterized by a scale. The first step is the pre-clustering of hadrons into 
hard final state jets and beam jets. One sets the hardness scale of the jets to d cnt . Then for 
every final state hadron hk and for every pair hk, hi one computes the corresponding value 
of the resolution variables dkB arid dki- There are several possibilities for the definition of 
the resolution variables. For instance, we may choose 

4b = v\, d\i = min(Pfe,Pj)i2L (4) 
where Rki is the distance in (y, 0)-space, 



Rki = - m) 2 + S - <Pi) 2 (5) 

and (pk, 6k, 4>k) are the cylindrical coordinates of the three-momentum of hadron hk with 
respect to the beam axis, with r]k = — lntan(6 , / c /2) being the corresponding pseudorapidity. 
(For other possibilities see ref. HU). Having calculated the resolution variables, one considers 
the smallest value among {dkB,dki}- If d^ is the smallest value and d^ < d cut , then hi and 
hj are combined into a single cluster with momentum according to a recombination 
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prescription and the algorithm starts again with hadrons hi and hj deleted from the final 
state and the 'pseudoparticle' of momentum added to the final state. If diB is the 
smallest value and diB < d cu t, then hadron hi is deleted from the final state and is included 
in the beam jets and the algorithm starts again. The algorithm stops if the smallest value 
is larger than the hardness scale d cut . At the end of the algorithm one has two beam jets 
and several hard final state jets. The second step of the algorithm is the resolution of the 
event structure into sub-jets. For this step one defines a resolution parameter y cut , y cnt = 
Qo/^cut — 1- Using this resolution parameter and the set of final state hadron momenta 
pre-clustered into the hard final state jets one performs a k± jet-clustering algorithm already 
familiar from studies in e + e~ annihilation. For the sake of simplicity, here we choose y cut = 1. 
With this choice we focus our attention to hard final state jets only, and suppress the second 
step of the clustering algorithm. 

The principle of parton-hadron duality implies that we use the same clustering algorithm 
at parton level as defined at hadron level. Thus the measurement functions used for the 
next-to-leading order perturbative calculation of N hard final state jet production are 

<SjV+l(Pl, • • • >PiV+l) = ®( d min 1} > d cut) + 9(^2^ < dcuOQ^S" 1 "^ > dcub) (6) 

and 

Sn(Pi, ■ ■ ■ ,Pn) = 9(42 > d cut ) (7) 

where 

d mL = mm ({P? ? dij} : i, j = 1, . . . , n, i=£ j), (8) 

and rf| 1 ^ n hl_ * Ar ' ) is the minimal value of the resolution variables after one clustering step. In 
eq. (|6|) the first term represents the iV + l-jet production, while the second one represents the 
production of N jets and either a soft parton, or a hard parton collinear with the beam axis, 
thus included in the beam jets or iV-jet production such that all final state partons are hard, 
but two of them are collinear, thus combined into a single jet. It is not difficult to check 
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that these measurement functions are infrared safe if any sensible recombination scheme 



TM is applied. In our analysis, we use the p t -weighted recombination. In this scheme the 



transverse momentum, pseudorapidity and azimuth of the pseudoparticle (ij) are defined as 

Pt(ij) =Pti+Ptj, (9) 

PuVti + PtjVtj nm 
vm = — , (10) 

, Pti<t>ti + PtjVtj 
<Pt(ij) = • (11) 

Ptm 

We remark that any other experimental cut, such as cut in the rapidity window, or a pt 
trigger should also be included in the measurement functions. In our analysis, we required 
that for the rapidity of jets \r)\ < 3 and the sum of the transverse momenta of the observed 
particles ^ um > 100 GeV. 

We now turn to the description of our results which were obtained at y/s = 1800 GeV 
machine energy and using the HMRS(B) [15[] parton distributions. In Fig. 1 we plot the total 
three-jet cross section both at Born level and at next-to-leading order for a fixed value of 
d C ut — 70 GeV vs which is the common renormalization and factorization scale. This plot 
demonstrates that over a wide range of fi values the theoretical uncertainty coming from the 
scale dependence is sizeably reduced in the next-to-leading order result as compared to the 
result of a Born calculation. In particular, one expects on general grounds that fi should be 
chosen of the order of the hard scale d cut . If one varies fi in the range of d cut /2 < fi < 2d CVLt , 
then the change in the next-to-leading order cross section is less than 1/6 of the change 
in the Born- level result. Thus the inclusion of the higher order correction decreases the 
uncertainty in the theoretical prediction by a factor of larger than 6. Similar conclusion can 
be drawn from plots at other d cut values in the range of 20 GeV < d cnt < 200 GeV. This 
can also be seen from Fig. 2, where the overall size of the three-jet cross section can be read 
off from a differential cross section G^ ut dcx/d<i cut plotted vs d cnt . The wide gray band shows 
the result of a Born level calculation with fi varied between d cnt /2 < fi < 2d cnt , while the 
narrower black band inside is the next-to-leading order result with the same scale variation. 

In order to give the reader some feeling about the errors of the Monte Carlo integrations, 



in Fig. 3 we plot the size of the Born-level cross section and the higher order correction to 
it together with the statistical errors of these results at fi = d cnt . The statistical error of the 
Born result is below 1 %, while the statistical error of the full next-to- leading order cross 
section plotted in Fig. 2 is below 10 %. 

In conclusion, we have calculated the three-jet cross section for the longitudinally invari- 
ant k± jet-clustering algorithm in hadron collisions for the simplified case when the matrix 
elements of all subprocesses are approximated by those of the pure gluon subprocess. The 
motivation for the particular choice of the jet definition is simply the pleasant property of 
the clustering algorithm that it uniquely assinges all final state particles to a certain jet. 
Using this definition one avoids the problem of jet separation in case of overlapping jets 
that occurs when cone jet definition is used. We used the subtraction method of ref. Jl2]| for 
canceling the infrared singularities. This method has the important feature that the physical 
quantity to be calculated is well separated from the theoretical problems of cancelation of 
infrared singularities and can easily be changed in a modular fashion in the Monte Carlo 
program. Thus the particular choice for the jet definition is by no means essential. We have 
found that the inclusion of the higher order correction improves dramatically our theoretical 
description of the three-jet cross section: the ambiguity coming from the arbitrary choice of 
the renormalization and factorization scales is reduced by a factor of at least six. Although 
the current analysis is not complete in the sense that we have not used the full QCD matrix 
elements, we anticipate similar conclusions once the effect of quarks is included. 
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FIGURES 
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FIG. 2. Differential three-jet cross section d^ ut da/dd cut vs d cut for 0.5d cut < \i < 2d cut at Born 
level (gray band) and at next-to- leading order (black band). 
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FIG. 3. Differential three-jet cross section d^ ut da/dd cut vs d cut for [i = d cut at Born level 
(crosses) and the higher order correction to it (histogram). The errobars indicate the statistical 
error of the Monte Carlo integration. 
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